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Anisotropic flow, which is sensitive to hydrodynamic expansion of the medium created in heavy 
ion collisions, has been extensively studied. Anisotropic flow measured using azimuthal correlations 
of flnal state particles is therefore contaminated by correlations unrelated to the reaction plane 
(nonflow). Currently anisotropic flow and nonflow cannot be separated experimentally in a model 
independent manner. Using PYTHIA simulations of p-hp collisions we show that nonflow approxi- 
mately factorizes in transverse momentum. This fact may be used to disentangle flow and nonflow in 
heavy ion data by performing a two-component factorization flt to Fourier coefficients of two-particle 
cumulants. 



I. INTRODUCTION 



The azimuthal distribution of particle momenta is an important observable in studying properties of nuclear matter 
created in relativistic heavy ion collisions. Pressure driven hydrodynamic expansion of the created medium transforms 
the initial configuration anisotropy of the colliding nuclei into the final state momentum anisotropy [T] . The initial 
spatial anisotropy decreases rapidly with the system expansion. Therefore the final state azimuthal anisotropy is 
primarily sensitive to the early stages of the system evolution i2, 'SI ■ 

Particle distributions are often quantified by the Fourier series in azimuthal angle (j) relative to the symmetry plane 

- — — ocl -I- V2t>„(pT)cos(n(0 -*„)). (1) 
dpTd© ■'^ 
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The symmetry plane is often called participant plane. It can be different from the reaction plane (which is defined 
by beam direction and impact parameter vector), and can be different for different n due to fluctuations Fourier 
coefficients w„ = (cos(n(0 — 4'n))) characterize the anisotropy for different harmonics. The second coefficient, V2, 
is called elliptic flow. It has been extensively studied and has provided compelling evidence for the creation of the 
strongly interacting quark gluon plasma at RHIC [SHE] . Studies of Vn of all harmonics are expected to yield valuable 
information about the properties of the quark gluon plasma [H ISHTT] . 

The participant plane is experimentally unknown and is estimated by the event plane ^ep which is constructed 
using all charged particles in an event [12] . Because of this, intrinsic particle correlations unrelated to the participant 
plane may contribute to the anisotropy. This contribution is generally called nonflow. Nonflow correlations include, 
for instance, jet correlations and resonance decays. 

Anisotropic flow can be also measured using the two particle correlation method. It was shown that the event 
plane and two-particle correlation methods are approximately equivalent . The two particle azimuthal correlation 
function can be expanded in a Fourier series in relative angle of the particle pair a,b: 

^3^Pair °o 

dp^HdA^ ocl + g ^VrMM cos(nA0) (2) 

In the case where correlations are only due to flow, particles are independently correlated to the participant plane, so 
VnA should factorize [Ti] : 

Ka(Pt,Pt) = vM)vnip'T) ■ (3) 

Thus Vn can be obtained for a given px by 

Vn'^HPT) = VVnAiPT,PT) • (4) 

Alternatively, if a reference Wn(PT'°^) known, w„ can be obtained by mixed-p^ pair: 
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In this paper we use w„ to stand for anisotropies measured by two-particle correlations (u„{2} as often used in 
literature [T5]). 

The flow factorization is a consequence of global correlations in an event: particles are all correlated to the common 
symmetry plane and particles are independent of each other. Nonflow, on the other hand, due to intrinsic correlations 
of particles, would not factorize a priori. Therefore two particle correlations in heavy ion events, including both flow 
and nonflow, would not generally factorize. However, a global fit to two-particle azimuthal correlations in Pb+Pb 
collisions using Eq. ^ yielded satisfactory results at low-pT- [IB]. This led to the suggestion that low-px correlations 
are dominated by anisotropic flow. 

In this paper we will show, by using PYTHIA simulations, that nonflow approximately factorizes, contrary to 
common perception. This is not surprising because jet fragments are all correlated to the jet axis, vE'jet, analogues 
to the symmetry plane in heavy-ion event. In the case of jet fragmentation, there is a conserved total momentum of 
thrust; analogously for the case of anisotropic flow, there is conserved total transverse momentum of zero. Except for 
this common correlation, jet fragments may otherwise be independent on each other. This would result in factorization 
of two-particle correlation in a jet: (cos(0i — (/)2)) — (cos(0i — 'I'jct)) (cos((/!'2 — ^jct))- This suggests that, although 
flow factorizes, factorization does not necessarily imply flow. In other words, the observed factorization in data alone 
does not necessarily lead to the conclusion that the correlations are dominated by flow. 

In the second part of the paper we will explore the possibility to extract flow and nonflow using a two-component 
factorization flt to mixed events of PYTHIA p+p collisions embedded in toy-model Monte Carlo simulations of 
hydrodynamical flow. 



II. NONFLOW FACTORIZATION 



One of the main sources of nonflow correlations in p+p collisions is jets. Particles in a jet are produced by 
fragmentation of a single high-pr parton, and are therefore correlated with the direction of the parton (jet axis). 
High-pT particles in a jet tend to be aligned with the jet axis while low-p^ particles have a lesser degree of coUinearity 
relative to the jet direction. Effectively, fragmentation process produces a "global" correlation with respect to the 
jet axis, and the correlation depends on pT- If jet fragments are indenpendent of one another, then nonflow (jet 
correlations) could also factorize. 

We investigate nonflow factorization using 210 M p+p events at y/s ~ 200 GeV generated using PYTHIA 6.4 [17] 
with default parameters except that we require a minimum momentum transfer of 3 GeV to ensure jet production in 
our simulations. Only stable charged particles are used in our study (pions, kaons, protons and antiprotons). The pt 
spectrum is shown in Fig. |3] (a) . 

In the experimental analysis of two-particle correlations an rj gap often is applied to reduce nonflow from resonance 
decays and intra (near-side) jet correlations. However inter (away-side) jet correlations are still present. We require 
A77 > 1 in our study: we correlate each particle from the backward pseudorapidity region (— 1 < 77 < —0.5) with 
particles at forward pseudorapidity (0.5 < 77 < 1). Two-particle correlation Fourier coefficients VnA are calculated 
using the direct Q-Cumulants method [TH1[T5], which avoids nested loops in the analysis. The Q- vector for a given 
Pt bin is defined by: 

N 
i=l 

where N is the particle multiplicity used in analysis. For two exclusive regions the two particle azimuthal moment 
IS given by Qn[PT)Q*n{PT)/{N{PT)N{PT)) whcrc iV(p^), N {p^j:) are multiphcities in those two T] regions. The two 
particle coefficient VnA is the event average: 

The error on y„A is simply a j \J N^vt , where a is sample standard deviation of VnA and A^ovt is number of events. We 
treat all errors as uncorrelated. 

Figure [ija) shows the two particle Fourier coefficient V2A obtained for PYTHIA p+p events. The correlation 
increases with px as expected as the particles are more closely aligned along the jet axis. Now we test applicability of 
the factorization hypothesis to correlation in p+p events. First, we calculate the single Fourier coefficient for nonflow, 
62, from the diagonal points of V2A distribution by Eq. Q. This is presented in Fig. [ijb). If factorization holds we 
expect that S2 obtained from V2A of mixed-pT pairs {pt,'Pt^) by Eq. ^ to be independent of p^°^ . Figure [ijc) shows 
single Fourier component (52 (pt) as a function of p^°^. In the case of lower px (1.5 — 2, 2 — 2.5), 82 is approximately 
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constant in the whole considered range of p^"^, while for the high pT (3 — 3.5 and 4 — 4.5 GeV/c) 52{pt) ~ const 
at pjP^ > 3 GeV/c. These results suggest that the nonflow correlations may approximately factorize in a limited pT 
range. 
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Figure 1: (Color online) (a) V2a(pt,Pt°') distribution in PYTHIA p+p. (b) 5^"'^ as a function of pr, (c) ^2 as a function of 
for several pr bins. 



In order to obtain a quantitative description of applicability of nonflow factorization, we fit V2a{Pt^Pt) data of 
PYTHIA p+p events with Eq. ([3|. Figure [2] shows the results. Because we expect approximate factorization to 
work best at high-pT, we cornpare fit results in different ranges of p^,p!^ pairs. We examine the fit quality for all 
Pt pairs (open circles in Fig. v2^, excluding pairs with the lowest pT only i.e. V2a{PttPt) points are excluded from 
fit if p^ < p^'" and p^ < p™'" (open triangles in Fig. and using only pairs with high-pT i-e. Px > P™™ and 



p^ > p™™ (open square and crosses in Fig. pi. The goodness of the fit is assessed with x^/NDF, where NDF is 
number of degree of freedom. In most cases the fit is driven by low-pT points and Eq. ^ underestimates V2A at 
high Pt- This may suggest that at low px there are additional processes that are not correlated with jet production. 
Fit quality (x^/NDF) improves when a higher threshold p™'" is applied, and the results are in good agreement with 
numerical data at high-pj- {pT > 3 GeV/c). Good fit quality at high-p^ is a convincing argument that nonflow due 
to jet fragmentation factorizes at higher px- 
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Figure 2: (Color online) (a) 52{pt) from Eq. ([sjl fit to PYTHIA p+p simulations in various pr ranges compared to 5^"^^ (see 
text for details), (b) Ratio of fit results to 5^'^^. 



III. NOVEL MEANS TO DISENTANGLE FLOW AND NONFLOW IN HEAVY ION COLLISIONS 



The relative contributions of flow and nonflow to the measured final state momentum anisotropy is an important 
question and several methods were proposed to distinguish them. Ollitrault et al. suggested that nonflow could 
be estimated using measurements in p+p collisions |20j . with the assumption that nonflow in heavy- ion collisions is 
identical to correlations in elementary p+p collisions. The PHOBOS experiment improved this approach by measuring 
azimuthal anisotropy using particles with a large rj gap, estimating long-range nonflow contribution using p+p data 
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and Monte Carlo event generators [5T]. Another approach was proposed in Ref. [H] where the measured correlation 
raw data were decomposed into elliptic flow and minijet correlations, where the minijet component was assumed to 
consist of a two-dimensional Gaussian in (A77, Ac/i) and a dipole moment in At/). In Ref. (TH], it was shown that if flow 
fluctuations are Gaussian, then it may be possible to disentangle flow, flow fluctuations and nonflow by simultaneous 
measurement of two-, four- and six-particle azimuthal cumulants. 

The approximate factorization of nonflow may provide us with a new means to separate flow and nonflow in heavy 
ion collisions and obtain their dependence. In a cluster model '23] where particles are composed of hydro particles 
and correlated particles from clusters, the final two-particle anisotropy for a given pT bin is 



V2a{PT,Pt) - [ -jfv2,hy + 2^ -^^^2,d I + 2^ ((cOS2A0)e;) . (8) 

cl / cl 



Here N^y is the number of hydro-medium particles, V2,hy - elliptic flow of hydro-medium particles, Nd - average 
number of clusters of one type, N^i - number of daughter particles in a cluster, J^d runs over all types of clusters, and 
V2,d is elliptic flow of cluster daughter particles which acquire flow due to anisotropy of clusters; (cos2A(/))c/ is the 
average cosine of twice pair opening angle in the cluster and is nonflow. Note that Eq[8] slightly differs from Eq. 22 
in Ref. [23] because we have assumed a Poisson distribution in number of clusters whereas in Ref. [23) the number of 
clusters was fixed. For a pair [p^,p^): Eq. (|8| can be generalized into 



where 



V2MyT) = V2{p^t)v2{p't) + E ^^^1^^ii^<^°^2A<^)^' (9) 



cl 



V2{Pt) = -^V2My + 2^ ^ V24 (10) 

cl 

represents the overall anisotropic fiow. If the nonflow term factorizes and only one type of clusters dominate, 

iV.iV.(p^)A^.(p^) ^ f0V^^(cos2(^^ rv/A^^(cos2(/-*. 

= S2{p^t)S2{p't) , (11) 
then V2A can be expressed as sum of two components, each of which factorizes: 

V2a{Pt,Pt) = MPt)MPt) + 62iPT)S2iPT) ■ (12) 



Here 82 is the "single particle" nonflow. 



/a^^(cos2 (0 , (13) 



where is a cluster axis. Note that in the literature [12] the symbol S is often used to stand for pair-wise nonflow, 
i.e. the second term of Eq. (l9| or Jl2|). In this paper we use 6 to stand for factorized "single particle" nonflow, and 
the pair- wise nonflow is S2{pt)S2{pj 



If nonflow factorizes, then Eq. ( 12 ) can be used in a two component fit to separate flow and nonflow. Note 
that the separation requires the pt dependence of flow and nonflow to be different, which should be the case in 
general. In order to study the possibility of separating flow and nonflow quantitatively, we consider a two-component 
model: hydrodynamical elliptic flow of bulk particles and nonflow (jet) correlations. The former component, which 
we refer to as "hydro", is approximated by a toy model simulation where we generate particles with V2 — a • 
where a = 0.1 GeV^^, which represents approximately the measured charged hadron elliptic flow for pT < 2 GeV/c 
in mid-central (30 - 40%) Au4-Au events at ^snn = 200 GeV jlj. Those particles are generated from soft pT 
spectrum dN/dpT oc pT-exp^^'^/^ where B ~ 0.23 GeV/c is obtained from fit to pion px spectrum at mid-rapidity 
reported by PHENIX in mid-central Au-|-Au events at y/s^N = 200 GeV (30-40% central events) [25]. For each event 
a random number of particles from the range of 46 - 133 is generated independently from a uniform distribution 
in two pseudorapidity regions: — 1 < 77 < —0.5 and 0.5 < 77 < 1. The multiplicity of generated Monte Carlo 
tracks corresponds approximately to the uncorrected multiplicity of mid-central Au-|-Au events at y^SNN = 200 GeV 
registered in the STAR Time Projection Chamber [26] , 
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Figure 3: (Color online) (a) Pt spectrum of flow ("hydro") and nonflow (PYTHIA p+p) event components, (b) u^'^^(pt) for 
pure hydrodynamical flow, PYTHIA p+p and mixture of both, (c) and (d): V2 = V2A(pT,Pr''')/w2'"'^(PT'^^) as a function of 
Pt°^ for pure hydrodynamical flow, PYTHIA p+p and mixture of both components for two bin of px of interest: 1.5 < pr < 2 
GeV/c (c) and 3 < pr < 3.5 GeV/c (d). 



We mix one PYTHIA event (nonflow correlation) with one "hydro" event. The pT spectra of hydro and PYTHIA 
particles are shown in Fig. [s] (a). The PYTHIA event is uncorrelated with the reaction plane of the hydro part of 
our events for this study. The "mixed" events are analyzed with the direct Q-cumulants method. 

Figure [s] (b) shows Fourier coefficient w^"'^ as a function of px calculated for pure "hydro" , PYTHIA and "mixed" 
events. The results are correlated with relative multiplicities of those components in an event: at low-p^, where 
"hydro" particles dominate, V2 follows closely the hydrodynamical flow while at high-px, where the jet component 
dominates, V2 is determined by nonflow correlations. Figure [s] (c-d) show V2 as a function of p^°^ for two regions 
of Pt- (c) 1.5 — 2 GeV/c, where hydro component dominates, (d) 3 — 3.5 GeV/c, where nonflow (jets) dominates. 
The observed dependence on reflects the relative contribution from flow and nonflow components: at low- 
Pt of interest, results are dominated by flow and V2 factorizes (is independent on p^^^) as expected. At higher 
p^^ V2 deviates from "flow-only" trend and then reassembles the nonflow correlations diluted slightly by hydro 
particles (see Fig. [3] (c)). However, V2 seems to factorize in the nonflow region as weU (Fig. |3] (c-d)). Note 
that in our mixed events the PYTHIA event is not correlated with the reaction plane of the hydro event. Adding 
in a correlation would change the details of the pjf'^ dependences in Fig. [s] (c-d). Moreover, even though the 
hydrodynamic flow dominates at low-px, it is slightly diluted by \ow-pT particles from the embedded p+p event. 
The contribution of hydro flow to the measured V2 scales with the ratio of number of hydro particles to the overall 
number of particles: V2{pt) = V2,hy{PT)Nhy{pT)/[Nhy{pT) + A^pythia(pt)], as in Eq. (10) (where W2,d = because 
our embedded PYTHIA event is random to the reaction plane). Similarly, the measured nonflow is given by (52 (pt) — 
S^^™^^{pT)NpYTmA{PT)/[Ni,y{pT) + A'^pythia(pt)] as in Eq. ^ where N^i = 1 in our case. 

We now fit formula ( |l2| ) to the V2A data of the mixed events. Fits were limited to pr < 5.5 GeV/c for "mixed" 
events to ensure sufficient statistical precision of the data. We expect that nonflow factorization works best at high-py. 
Therefore we compared flts in three different cases: (a) < < 5.5 GeV/c and < < 5.5 GeV/c, (b) 2 < < 5.5 
GeV/c or 2 < < 5.5 GeV/c i.e. points with low px pairs are excluded from fit, (c) 2 < p^ < 5.5 GeV/c and 
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Figure 4: (Color online) Fitted flow «2 and nonflow S2 as a function of pr compared to the expected input flow and nonflow 
(see text for details). 



2 < < 5.5 GeV/c i.e. only high-pT points are used in the fit. In each case fit quality is good: /NDF sa 1. 
Goodness of the fit at \aw-pT for mixed events is better than in the case oi p+p data alone because low-p^ is dominated 
by hydro flow in the mixed events, which factorizes. In order to compare the fit resuhs with expected values of V2 and 
82, we present the ratio of V2 and (S2, scaled by relative muhiplicity, to the input parameters (Fig. |4](c) and (d)). In 
the case of non-fiow, there is a 30 — 45% deviation at lowest-pr; at pt > 2 GeV/c good agreement is observed. This 
is expected because we showed in Sec. |ll]that nonflow factorization works best at high-p^-. The fit result for flow 
is also good - there is a small deviation (less than a few %) in the region where data provides sufficiency statistical 
precision {pT < 3 GeV/c). 

We note that the fit of a quadratic function of the form of Eq. ( 12 ) to the V2A data is a mathematical problem with 
possibly multiple solutions. Different solutions can result from a different set of initial values of the fit parameters. 
However the solution representing the real physical correlations can be distinguished from the purely mathematical 
ones when the high px region is considered, where the nonflow correlations dominate and anisotropic flow is expected 
to be small. 

It is worthwhile to examine quantitatively the higher harmonic anisotropic flow and two-particle azimuthal corre- 
lations measured in Pb-|-Pb collisions [Tnill71[5S|. In Ref. 157] the higher harmonic t;„{2} (n = 2,3,4,5) were first 
obtained from two particle correlations, where particles of interest a with a given p^ were correlated with reference 
particles c of all pt- The azimuthal correlations between triggered and associated particles were then compared to cor- 
relation shapes expected from the measured f2{2}, t'3{2}, W4{2}, W5{2} harmonics evaluated at the pr corresponding 
to trigger and associate particles. It was shown that the measured Fourier coefficients describe the trigger-associate 
correlations. Let us examine what difference would be expected between the trigger-associate correlations and the 
composition of products of the Fourier harmonics measured with a third particle c. In the case of two-component 
factorization, Vn{2} for the trigger particle 6 is given by v'^ip^) = (t;„(p^)t;„(p^) -l-<5„(py)5„(py))/i/v2(p^) -|- S^{p^), 
where the denominator is the c particle harmonic measured by two-particle correlation with both particles from all 
Pt- Likewise for the anisotropic flow v'^^{p^) of associate particle a. The deviation of trigger-associate correlation 
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from composition of single particle harmonics would be 



Vn[PT)Vn{PT) + 5n{PT)^n{PT) 



1 



(14) 



in case of small nonflow relative to flow. Another approach is to take the trigger (and associate) particle harmonic 
directly from two-particle correlations in the same trigger (associate) px region |28| . The deviation in this approach 
would be 



Ka(p^,P^) 



- 1 = 



Vr^Aip'^T^PT) 



vMXip'r) ^V,M,p^T)^/VnAip'T,p'T) 



1 



2 KVniPr) Vn{PT) 



(15) 



As seen from Eqs. (14| and (15l, even for relatively large nonflow Snipr) /vnipr) ~ 10%, the deviation amounts to 
only the order of 10"'^ and would be very challenging to observe with current statistical precision of the data, which 
is roughly what the results in Refs. [571 [IH] indicate. Therefore it is not surprising that global fit in Ref. [T^ gave 
satisfactory results at low-p^- If ^n{PT) / "^niPT) — 5n{p^) / Vn{p\') i-e. flow and nonflow have the same pt dependence, 
then V^nA (Pt,Pt) factorizes precisely even though nonflow is present. 



IV. SUMMARY 



We have shown by PYTHIA studies of p+p collisions that nonflow correlations due to jet fragmentation approx- 
imately factorize. Based on this observation we developed a novel method of fitting a two-component factorization 
expression to heavy ion data in order to extract anisotropic flow and nonflow. We used PYTHIA events embedded in 
toy-model Monte Carlo simulations of hydrodynamical flow to demonstrate that this method works well for a mixture 
of flow and nonflow. 

We focused on the second harmonic (y2 and ^2) in this paper. In principle, a two-component factorization fit 
should work for other harmonics as well, although the sign of nonflow contribution has to be established for specific 
case. For instance, jet correlation from PYTHIA would give a negative contribution and therefore V3A (Pti-Pt) = 
v?,,hy{pT)'"i.hviPT) - S3{Pt)^3{Pt) should be used in the fit. 
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